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THERMODYNAMICS  OF  SURFACE  CLUSTERS  - DIRECT  OBSERVATION  OF  Re.,  ON  *(211)' 

by 

Kaj  5coLc('3\  John  D.  Wrigley,  and  Gert  Ehrlich 
(c) 

Coordinated  Science  Laboratory  , Department  or  Metallurgy,  and  Department 

of  Physics 

University  of  Illinois  at  Urbana-Champaign 


ABSTRACT 

Direct  observations  in  the  field  ion  microscope  make  it  possible  to 
determine  the  ratio  of  bound  to  dissociated  dimers  at  equilibrium  on  a 
crystal  surface.  In  principle,  the  thermodynamics  of  dissociation  are 
accessible  from  the  temperature  dependence  of  this  ratio.  To  establish 
the  feasibility  of  such  measurements,  the  kinetics  governing  changes  in  the 
dimer  probabilities  are  worked  out.  This  is  done  for  isothermal  processes, 
and  for  conditions  typical  of  the  quench  from  the  temperature  at  which  the 
equilibrium  distribution  is  established  to  the  temperature  at  which  observa- 
tions are  made  in  the  field  ion  microscope.  Techniques  are  developed 
to  correct  tor  distortions  of  the  equilibrium  distribution  during  the 
quench.  These  rely  upon  measurements  of  the  ratio  of  dimers  in  configura- 
tion 1 and  0 at  low  temperatures,  combined  with  observations  of  ? , the 
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fraccion  dissociated  at  high  temperatures.  Experiments  conforming  to  this  I 

scheme  have  been  carried  out  for  rhenium  dimers  on  W(211).  They  yield 

internal  energies  of  -3.7  +1.1  kcal/mole  for  rhenium  dimers  in 

configuration  1,  and  -2.3  + 1.2  kcal/mole  for  dimers  in  configuration  0, 

measured  relative  to  the  energy  of  dissociated  pairs.  Our  results  suggest 

important  contributions  from  indirect  interactions. 


3 


During  Che  past  decade,  direct  observations  in  the  field  ion  microscope 
have  yielded  considerable  quantitative  information  about  the  behavior  of 
individual  metal  atoms  diffusing  on  crystal  surfaces.1.  More  recently, 

2 

these  observations  have  been  extended  to  clusters  of  interacting  adatoms, 

and  for  dimers  we  now  have  available  a fairly  detailed  description  of 

3-5 

surface  transport  on  the  acoraic  level.  The  emphasis  of  quanti- 
tative investigations  has  so  far  been  upon  the  kinetics  of  atomic  events. 

In  contrast,  the  energetics  governing  the  association  of  metal  adatoms 

into  clusters  have  not  been  extensively  explored.  A theory  of  the  forces 

6 -8 

between  adatoms  on  metal  surfaces  is  beginning  to  emerge,  and  some 

9 

observations  of  cluster  dissociation  have  been  reported.  However,  we  do 
not  yet  have  a quantitative  piccure  of  surface  interactions  between  metal 
atoms,  although  these  are  fundamental  to  an  understanding  of  cohesion  in 
crystal  layers. 

In  this  paper,  we  address  ourselves  to  the  first  quantitative 

determination  of  the  thermodynamics  governing  dissociation  of  dimers 

on  a crystal  surface.  The  primary  emphasis  is  upon  the  conditions  and 

procedures  appropriate  for  such  thermodynamic  studies.  We  concentrate 

upon  one  illustrative  example  - rhenium  on  the  (211)  plane  of  tungsten  - 

to  provide  an  indication  of  atomic  forces  in  surface  clusters.  Much  is 

4 

already  known  about  the  chosen  system  from  previous  studies.  As  appears 
from  Fig.  1,  the  (211)  plane  is  a channeled  surface  on  which  atomic  motion 
is  confined  to  one  dimension.  This  simplifies  both  the  experimental 


observations  and  their  analysis.  When  two  atoms  are  placed  in  adjacent 


in'!  plane  of  tungsten.  Dime 
hannels  cf  close-packed  atoms 
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channels  of  this  plane  and  equilibrated,  the  probability  of  finding  them 
dissociated  rather  than  bound  to  each  other  provides  a measure  of  the  free 
energy  change  on  dissociation;  studies  at  different  temperatures  should  then 
yield  both  the  energy  and  entropy  of  dissociation. 

Experimental  determinations  of  this  probability  ratio  are  quite 

10  _ . - 

straigntrorvara  in  tae  tieic  ion  microscope.  me  sample  surrace  is 

kept  at  the  desired  temperature  for  a time  interval  sufficient  to 
establish  thermal  equilibrium,  and  then  quenched  to  cryogenic  temperatures. 

The  distribution  of  adatoms  over  the  surface  is  observed  at  low  temperatures 
only,  to  avoid  any  possible  perturbation  of  the  system  by  the  act  of  observa- 
tion. It  is  assumed  chat  the  observations  provide  a close  approxima- 

» 1 

tion  to  Che  equilibrium  distribution  at  the  chosen  temperature.  This 
assumption  is  critical.  If  the  cooling  rate  is  not  high  enough,  extensive 
rearrangement  of  Che  adatoms  may  occur  during  Che  quench,  invalidating 
the  measurement.  In  fact,  attaining  high  quenching  rates  in  the  field 
ion  microscope  presents  quite  3 serious  experimental  problem.  Even 
though  our  main  concern  is  with  equilibrium  properties,  we  must  therefore 
also  analyze  the  kinetics  governing  the  redistribution  of  dissociating 
dimers . 

After  a brief  presentation  of  the  thermodynamic  formalism  for 
dissociation  and  of  some  experimental  particulars,  we  examine  rate  processes 
tor  dissociating  dimers  in  Sec.  III.  The  kinetics  governing  changes  in 
the  probability  distribution  of  dimers  on  a surface  are  worked  out  under 
isothermal  conditions,  and  are  then  adapced  to  treat  changes  during 
quenching.  Based  on  this  analysis,  procedures  are  outlined  (in  III-C)  to 
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3llow  tor  the  distortion  of  the  high  temperature  distribution  on  cooling. 
These  procedures  make  possible  meaningful  studies  of  the  thermodynamics 
of  dimers  on  a surface.  Actual  thermodynamic  measurements  on  rhenium 
dimers,  which  constitute  the  central  theme  of  this  effort,  are  presented 
in  Sec.  IV  and  V. 

I.  DISSOCIATION  EQUILIBRIUM  OF  DIMERS 
At  the  temperatures  considered  here,  adatoms  on  the  (2111  plane  of 
tungsten  are  confined  to  channels  formed  bv  surface  atoms  close-packed  in 
the  .111’  direction.  When  two  adatoms  are  placed  in  adjacent  rows  of 
such  a surface  and  equilibrated,  Che  probability  of  finding  the  atoms 
at  a distance  i from  each  other  ls^ 


Pi  " C(2-$l0)(L-l)exp  - Wt/kT. 


(1) 


Here  _i  denotes  the  separation  between  adatoms,  measured  in  the  direction  of 

the  channels,  in  units  of  the  spacing  i along  the  .111];  L is  the  number  of 

sites  accessible  to  an  atom  in  a channel,  5.,  is  the  Kronecker  delta,  and 

lO 

C is  a norma licat ion  constant,  to  insure  that 


L-l 
£ P. 
i-0  L 


(2) 


Interactions  between  the  atoms  enter  through  W^,  the  potential  of  mean 
force;  in  our  situation,  it  can  be  equated  to  the  change  in  the 
Helmholtz  free  energy  F on  moving  two  adatoms,  initially  at  an  infinite 
separation,  to  a distance  _i  from  each  ocher.  Since  the  normalization 
constant  C itself  depends  upon  the  temperature,  measurements  of  ? as  a 


function  of  the  temperature  T do  not  suffice  to  determine  the  interaction. 
Instead,  it  is  necessary  to  cake  the  ratio  c two  probabilities  - one  at 
distances  so  Large  that  interactions  between  adatoms  are  negligible,  the 
other  at  the  separation  JL  of  interest.  Based  on  our  previous  work,  we 
will  assume  that  once  the  distance  between  adatoms  equals  or  exceeds  two 
s pacings  along  the  'til-  » interactions  have  dropped  to  zero  and  the  two 
adatoms  are  free  (W,^„  = 0).  Anv  oair  for  which  t < 1 is  counted  as  bound. 

The  probability  PQ  that  a pair  be  dissociated  is  now  given  by 

L-l 

P - E P.  - (L-l)(L-2)/2,  (3) 

U i-2  1 

and  the  ratio  P^/P^  is  just 

P./PD  = [2/ (L-2)jexp  - W^kT 

“ [2/(L-2)]exp  S^/k  exp  - E,/kT  . (4) 

Here  and  throughout  this  paper  all  thermodynamic  quantities  are  referenced 
to  the  properties  of  dissociated  pairs;  the  incernal  energy  and  entropy 
of  dissociated  dimers  3re  therefore  set  equal  to  zero.  A semi- 
logarithmic  plot  of  Pj_/?D  versus  1/T  snould  yield  both  the  internal 
energy  and  the  entropy  of  pairs  in  state  1,  provided  the  number  of  sites  L 
in  the  diffusion  channel  is  known.  It  must  be  noted,  however,  that  apart 
from  the  statistical  problems  usual  in  quantitative  studies  of  individual 
atoms,  there  is  now  the  additional  difficulty  of  establishing  the  temperature  T 
to  which  the  observed  probabilities  pertain.  This  will  be  discussed  in 


Sec.  III. 
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II.  EXPERIMENTAL  CONDITIONS  - QUENCHING  RATES 
Our  experiments  hinge  on  the  ability  to  distinguish  dissociated  from 
bound  dimers.  As  indicated  in  Fig.  2,  this  is  easily  done  in  the 
field  ion  microscope.  The  procedures  and  equipment  used  to  obtain  field 
ion  images,  and  especially  to  establish  the  separation  between  atoms  in  a 
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recording  uf  images 


has  been  improved,  however,  by  using  an  Olympus  0M2  camera,  with  f/1.2 
optics,  which  automatically  adjusts  exposures  to  the  correct  time  setting 
for  Daylight  Ektacnrome  film. 

Most  important  for  the  present  work  is  the  extent  to  which  the 

distribution  of  bound  and  dissociated  pairs  is  altered  by  quenching  to 

the  cryogenic  temperatures  suitable  for  field  ion  microscopy.  Such 

rearrangement  depends  upon  how  rapidly  the  sample  cools  after  equilibration. 

This  temperature  decay  has  been  measured  for  the  two  support  loops  used  in 

this  study.  Both  are  made  of  .007  in.  diameter  tungsten;  the  end-to-end 

12 

length  is  *»  3 cm  for  the  first,  and  only  =>  2 ca  for  the  second.  During 
the  quench,  the  resistance  of  the  support  is  monitored  by  a 4-lead 
technique.  The  emf  across  the  potential  leads  to  the  loop  is  displayed  on 
an  oscilloscope,  as  in  Fig.  3.  It  is  thus  possible  to  follow  rapid  changes 
in  the  temperature  of  the  support,  which  may  exceed  one  hundred  degrees 
per  second.  The  time  for  the  support  to  cool  down  to  300°K  is  critical 
for  freezing  in  place  the  high  cemperature  distribution.  During  this  interval, 
the  cooling  curve,  shown  in  Fig.  3,  can  be  approximated  by  the  hyperbolic 
re lation 


1/T  = 1/T„  + bt 


(5) 


1 


Configurations  of  bound  and  dissociated  rhenium  dimers  on  W(211) 
Hard-sphere  model  on  left,  field  ion  image  on  right.  Dimer 
configurations,  from  top  to  bottom:  0,  1,  2,  and  4:  the  top  two 
are  bound,  the  bottom  two  free.  Interatomic  separation  of  atoms 
in  dimer  established  by  techniques  in  reference  4. 
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The  decay  race  parameter  b is  essentially  independent  of  T , the  initial 

-4 

temperature  ot  the  support  loop.  For  the  long  support,  b » 5.75  x 10 

-1  *3  -1 

(deg  sec)  ; for  the  shorter  loop,  b - 3.0  x 10  (deg  sec) 

With  the  cooling  race  established,  it  is  now  possible  to  assess  the  extent 

to  which  the  original  distribution,  established  in  conformity  with  Eq.  (4) 

by  equilibration  at  a high  temperature,  is  affected  during  the  quench  to 

the  low  temperatures  at  which  observations  are  made. 

III.  KINETICS  OF  DIMER  FORMATION  AND  DISSOCIATION 
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To  understand  how  an  equilibrium  distribution  of  dimers  is  perturbed  by 
changes  in  the  temperature,  the  formalism  for  the  kinetics  of  dimer 
dissociation  at  constant  temperature  is  examined  first.  In  III.B  we  then 
analyze  the  effects  of  quenching  on  an  equilibrium  distribution  of  model 
dimers.  Finally,  the  implications  of  this  kinetic  analysis  for  studies  of 
Che  dissociation  equilibrium  are  examined,  and  techniques  are  devised  for 
utilizing  experimental  observations  in  such  a way  as  to  circumvent  changes 
in  the  actual  high  temperature  distributions  during  cooling. 

A.  Isothermal  Dissociation 
1.  Formalism 

Migration  of  dimers  in  an  infinite  one-dimensional  channel  is  readily 

described  by  noting  the  position  of  the  center  of  mass,  together  with  the 
, 3 5 

configuration  of  the  dimer  ’ ; th*»  latter  is  given  by  the  separation  _i 
between  Che  two  adatoms  in  che  pair  (measured  in  units  of  the  atomic  spacing 
l along  the  channel).  The  conventions  governing  the  rate  constants  for 


wmm 
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acorn  jumps  from  one  site  Co  a specified  neighboring  site  are  indicated 
schematically  in  Fig.  4.  The  rate  is  denoted  by  a if  the  jump  increases 
the  pair  separation,  and  by  9 if  it  causes  a decrease.  The  subscript  to 
the  symbol  indicates  the  configuration  prior  to  the  jump;  jumps  can  occur 
only  between  adjacent  sites. 


The  probability  p that  the  center  of  mass  of  a dimer  in  configuration 

X,  L 
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i^  be  at  position  x therefore  changes  in  accord  with 


dp 


xj. 


i(P„ 


dt  »x,i  1-1  x-1 , i-1 


+ Px+l.i-l*  “ 2(3  i + Qri)px,i 


+ 3i+l(Px-l,i+l  + px+l , i+1 ^ 


x =*  0,  +1 , +2,  ...  , (6) 


where  =*  3q  =■  0.  We  are  not  really  concerned  with  the  position  dependence 


of  the  probability,  however;  our  emphasis  is  rather  on  the  probability  P. 


of  finding  a dimer  in  a given  state  i_.  By  definition, 


P . = E p 
l x *x,i 


(7) 


Summing  Eq.  (6)  over  all  values  of  x now  gives  the  time  dependence  of  the 
probability  P^: 

*1  " 2^i-lPi-l  “ (afi  +3i)Pi  + 3i+lPi+lJ 

i = 0,1,2,... 


(8) 


In  a finite  system  the  number  of  accessible  configurations  is  limited. 
Henceforth, we  confine  ourselves  to  a set  of  K configurations,  ranging  from 


■ — I 


Fig.  4.  Race  constants  for  jumps  of  dimers  in  one-dimensional  diffusion. 


Center  of  mass  positions  are  indicated  by  crosses  and  labelled 
with  italic  numbers.  Arrows  show  atomic  jumps  allowed  between 
nearest  neighbor  positions.  Jump  rates  to  higher  configuration 
are  given  by  y,  to  lower  configuration  by  3;  subscript  labels 
the  starting  coni iguration. 


i ■ 0 :o  i » K-l,  and  assume  that  there  is  no  interchange  with  higher 
configurations,  ao  that  a,,  , * 5 „ - 0.  With  this  assumption,  we  are  Lett 
with  the  equations 


? ■ '»"  -v  ? ^ a i 

3 0 0 ‘l*  l' 


O m P 
• 1 “i 


2fjv  P - (.y  > P +.  i o ' 

1 0 0 '‘l  'Ll  2*2* 


PJ  " :Corj-l?J-l  ' * V?J 

PK-l  ' ^K-2?K-2  • 3H?H:  * 


J - -.3 K-2 


The  solution  Co  chis  sac  of  differencial  relations  is  well  known;  it 
corresponds  Co  the  concentration  of  the  iCh  incarmediace  found  in  a aaquanca 
of  X coupled,  first-order  chemical  reactions,  and  is  given  by 


?i  - ,E,  Vu  - V • 

X«l 

The  \.  are  eigenvalues  of  the  matrix  2>t,  where  M is  defined  as 


-orQ  (orL  + 3 L) 


(or,  + 3-) 


«<V  ( y ♦ 3 N — ^ 

K-2  k--  X-l 


-a..  , 

K-  - 


2he  coefficients  A.  are  given  bv  the  comoonenta  of  the  correspond ing 

LX  3 


eigenvector . 


3och  the  eigenvalues  and  eigenvectors  can  be  evaluated  by 
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However,  for  a channel  of  L sites  we  know  from  Eq.  (1)  that  when  i ^ f, 
and  hence  VT  * 0, 

Pi+1/Pl  " i»£. 


lb 


To  satisfy  this  requirement,  the  probabilities  given  by  Eq . (10)  are 
weighted  by  a).,  where 


x.  « (L-i)/  I (L-i)P  . 

i*0  L 


Reed  has  demonstrated  by  Monte  Carlo  simulation  that  this  procedure  provides 
a satisfactory  approximation  to  dimer  dissociation  on  finite  channels,  even 
far  removed  from  the  steady  state.  The  solutions  given  by  Eqs.  (10), 
weighted  in  accord  with  Eq.  (13),  will  henceforth  be  used  to  describe 
the  probabilities  of  model  dimers  at  constant  temperature. 


2.  Time  Evolution  of  Probabilities  for  Model  Dimers 

In  order  to  show  how  the  probabilities  P^  approach  their  steady  state 

values,  Eq.  (10)  has  been  evaluated  for  a model  system,  described  in  Fig.  5. 

This  system  is  patterned  loosely  on  the  properties  expected  for  Re,,  on 

the  (211)  plane  of  tungsten,  as  suggested  by  previous  studies.  In  the  model, 

interaction  effects  are  neglected  beyond  configuration  1;  that  is,  dimers 

with  an  interatomic  separation  of  2 spacings  are  assumed  to  behave  as  free 

acorns.  The  jump  rate  o'  of  free  atoms,  and  the  rate  3 ^ at  which  dimers 

jump  into  the  0 configuration,  are  sec  equal  to  the  values  found  in 
4 

experiments  with  rhenium  on  W(211).  The  temporal  evolution  of  the 
probabilities,  shown  in  Fig.  6,  is  specific  for  the  model  system.  However, 
the  particular  choice  of  the  race  constant  3^  does  not  significantly 
affect  the  outcome;  a change  in  3 by  a factor  of  two  does  not  perceptibly 
change  the  plots,  provided  the  ratio  o^/3 ^ remains  fixed. 


Interatomic  ceoarat'on  > J! 


Jimer  Configuration 


'ig.  5.  Schematic  potential  diagram  for  model  system.  Dimers  in  configura- 
tion 2 and  higher  are  completely  dissociated;  their  potential 
energy  is  taken  as  the  tero  of  reference.  Rate  constants  are 


assigned  the  following  values:  a./3.  - 0.4  exp  1234.o,'R,T; 
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x 
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or  “ 2.95  x IQ*“  exp-19 ,300/R,T  sec 

u 

experimental  guantities  for  rhenium  on  V(211) 
R„  ■ L.987  cal  (deg  mole'  " . 


. - 3 x v exp-3500/R,T,  3 » 2.5  x L0*“ex?-13 ,300/R,T  see"1, 

in  u l u 

he  last  two  rates  are 


The  t tine  dependence  of  the  probabilities,  which  Is  displayed  in 


Fig.  o for  a dimer  starting  in  the  0 configuration,  is  complicated  and 
has  contr ibutlons  from  a number  of  exponential  terms.  However,  equilibra- 
tion between  states  0 and  1 occurs  instantly  on  the  time  scale  appropriate 
to  dissoc lation.  Only  .02  seconds  after  starting  in  the  0 conf igurat'on, 


the  ratio  is  within  0.1*1  of  the  equilibrium  value;  a plot  for  dimers 

starting  in  the  l conf iguration  is  indistinguishable  from  the  one  shown.  The 
equilibration  of  the  higher  configurations  takes  considerably  longer,  and 
follows  a course  typified  by  the  probability  ?.,  - a rapid  initial  rise,  in 
which  the  equilibrium  value  is  overshot,  followed  by  a gradual  diminution, 
as  redistribution  to  the  ocher  conf igurations  occurs.  The  sum  total  of 
che  probabil  it  ies  for  states  other  than  0,  l or  2 is  included  under  the 
symbol  ; this  quantity  (and  also  the  probability  PQ  of  finding  a dimer 
dissociated^  rises  slowly  but  monotonical ly  to  its  long  term  value. 

3.  kinetics  of  Quenching 

Of  greater  interest  for  this  study  are  changes  in  the  distribution  which 
occur  after  equilibration,  while  the  surface  cools  to  the  temperature 
at  which  observations  are  made.  These  changes  are  evaluated  by  approxi- 
mating che  actual  hyperbolic  cooling  curve  of  the  sample  with  a series  of  1 

isothermals,  A degrees  apart,  each  terminated  by  an  ins tantaneous  drop  to 
the  next  lower  temperature.  This  procedure,  suggested  schematically  in 
Fig,  7,  makes  it  possible  to  determine  tne  extent  of  distortion  during  the 
quench  without  extending  the  formalism  of  the  previous  section. 

Estimates  of  such  distortions  are  shown  in  Fig.  3 for  our  model  system 
after  equilibration  at  four  temperatures  at  which  significant  dissociation 
occurs.  From  these  curves  it  is  evident  that  che  procedure  outlined  in 

1 i 


Fig.  7.  Representation  of  cooling  curve  for  sample  support  by  isothermal 


intervals  separated  by  step  function  drop  of  ! degrees.  For 
quantitative  estimates  of  distortion  during  cooling,  2-4° 


I 


Simply  observing  Che  ratio  ? /?  after  equilibration  at  different  tempera- 
tures will  not  work:  the  distribution  obtained  by  equilibration  at  high 


temperatures  is  significantly  altered  during  the  quench  to  the  temperature 
of  observation,  even  for  equilibration  at  temperatures  as  low  as 
350° K.  In  large  measure  the  problem  arises  from  the  rapid  interconversion 
between  states  0 and  1.  As  is  apparent  in  Fig.  9,  the  conversion  is 
significant  even  at  temperatures  below  350°K;  only  below  300°K  is  the  ratio 

?,/P.  frozen  in. 

1 0 

The  probability  of  finding  a dimer  dissociated  is  much  better 
behaved.  Dissociation  is  the  limiting  step  in  establishing  an 
equilibrium  distribution,  and  stops  quickly  during  a quench.  The 
reverse  process,  recombination  into  a bound  pair,  faces  two  obstacles: 
che  comparatively  high  barrier  to  diffusion  for  individual  atoms  in 
this  system,  and  the  need  on  the  average  to  diffuse  ever  an  appreciable 
distance  before  encountering  Che  other  adatom  in  the  adjacent  channel. 
During  a quench,  these  effects  preserve  the  number  of  dimers  that  had 
dissociated  at  high  temperatures.  Sven  after  equilibration  at  425°K,  the 
probability  p is  changed  by  only  =»  12.5%  during  the  subsequent  cool-down. 
To  the  extent  that  we  focus  upon  this  quantity,  valid  thermodynamic 
information  should  therefore  be  attainable. 

C.  Analysis  of  Dissociation  Data 
3ased  on  these  calculations  for  model  dimers,  it  is  now  possible 
to  devise  a scheme  for  properly  analyzing  real  systems.  In  order  to 
derive  the  thermodynamics  of  dissociation  from  Eq.  (4),  we  require 
reliable  data  for  ?^/P  . 


Only  ?D , however,  can  be  derived  from 
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measurements  at  high  equilibration  temperatures  without  major 
distortion  during  quenching.  We  therefore  estimate  values  of  P^  for  the 
temperatures  at  which  dissociation  is  significant  by  extrapolating  P^/P^ 
from  lower  temperatures,  at  which  distortion  of  this  ratio  is  not  too 
severe.  These  estimates  are  straightforward,  provided  sound  data  are 
available  for  the  thermodynamics  of  dimers  in  state  1 relative  to  state  0. 

Consider  an  experiment  on  a channel  of  length  L,  for  which  a total  of  M 
observations  is  available  at  a given  equilibration  temperature.  Dimers 
are  in  state  l in  of  these  M observations,  they  are  in  state  0 in  Mg 
observations  and  are  dissociated  in  M^.  Only  M and  give  significant 
information  about  dimer  behavior  at  the  equilibration  temperature  : the 

C 

population  of  states  0 and  1 changes  coo  much  during  quenching.  Never- 
theless, we  can  get  Mg  and  M^  indirectly,  by  noting  that 

Mq  - iM-Mgj/Cl+f^) 

M1  - (M-MgJR  /(1-tt^  . (14) 


Here  the  ratio  M^/Mg  for  a channel  of  L sites  is  denoted  by  R^.  As 
already  indicated,  this  ratio  is  derived  by  extrapolation  from  low 
temperatures.  Thus,  from  Eq.  (14)  we  can  estimate  Mg  and  M for  the  high 
temperatures  at  which  dissociation  is  favored;  this  can  be  accomplished 
without  the  distortions  inevitable  in  an  attempt  to  observe  the  distribution 
after  quenching  to  cryogenic  temperatures. 


The  program  for  determining  the  thermodynamics  of  dimer  dissociation 
now  clear.  It  requires  measurements  of  Pg,  the  probability  of  finding 


1 
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dimers  dissociated  after  equilibration  at  high  temperatures  (in  the 
o 

vicinity  of®’  400  K).  Also  needed  are  accurate  values  of  P^/P^  for 
lower  equilibration  temperatures  (300°K  and  less),  to  allow  extrapolation 
of  this  ratio  to  temperatures  at  which  dissociation  is  favored. 


I 


IV.  EQUILIBRIUM  BETWEEN  RHENIUM  DIMERS  IN  STATES  0 AND  1 

Measureraep.es  of  Che  relative  population  of  rhenium  dimers  present  in 

4 

states  1 and  0 have  previously  been  reported  for  W(211).  These  studies 
were  carried  out  over  a wide  range  of  temperatures,  extending  up  to  390°K. 


For  the  higher  temperatures,  however,  considerable  distortion  of  the 
distribution  must  have  occurred  while  the  sample  cooled  to  the  temperature 


of  observation.  In  fact,  above  330  K,  the  previous  data  on  P^/PQ  (in  Fig.  14 
of  reference  4)  appear  insensitive  to  the  equilibration  temperature;  this  strongly 
suggests  that  annealing  during  the  quench,  rather  than  the  equilibration, 
dictated  the  observed  distribution.  More  reliable  values  of  R = P^/P^ 
are  therefore  needed  to  allow  extrapolation  of  this  ratio  into  the 


dissociation  regime. 


A.  Observations 


Additional  measurements  of  rhenium  dimers  in  states  0 and  1 have  been 
made  at  temperatures  from  270°  to  327°K.  The  observations  are  summarized 
in  Table  I,  and  are  plotted  in  Fig.  10.  Only  statistical  errors  are 


indicated,  derived  from  the  standard  relation' 


var  R = R-/[(M-l)PnP.] 


this  assumes  that  dissociation  is  negligible  at  these  low  temperatures. 

On  a channel  of  L sites,  R^,  the  ratio  of  pairs  in  state  1 to  state  0, 
is  related  to  the  thermodynamic  properties  of  the  two  states  by 


\ ' (Pl/P0)L  M 2^L-1'»/I']exp(S1-S0)/k  exp  - (E^E^/kT  . 
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Fig.  10.  Temperature  dependence  for  ratio  R of  rhenium  dimers  in  state  1 
compared  to  state  0,  on  W(21L).  Raw  data  on  bottom  (ordinate  a 
left)  yields  thermodynamic  daca  in  sat  (A).  Values  on  top 
corrected  for  distortion  during  quench  (ordinate  at  rights,  giv 
daca  sec  (3).  Error  bars  indicate  standard  deviation  arising 
from  statistical  affects;  curves  obtained  by  least-squaras  fit 
to  daca  poincs  weighted  by  R'.'var  R. 
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or  'll'  :.:x 


In  che  experiments,  dimers  are  iniciaLLy  deposited  in  che  cancer  of  rhe 

channeLs.  Ac  che  Low  cemperacures  ac  which  states  0 and  L are  equili- 

braced,  dimers  rarely  move  co  che  edge  of  che  plane;  whan  chay  do,  che 

observacions  are  disregarded.  Adjustments  for  che  Lengch  L of  che  channels 

are  cherefora  noc  made  in  Eq.  (16)  - che  racio  M.  /M.  derived  from  che  measure- 

1 u 

nencs  is  creaced  as  equivalent  co  Ra,  che  value  of  ? /?^  on  an  infinice 
channel.  A semi-logar ichtnic  ploc  of  che  observed  values  of  ? /P^  againsc 
reciprocal  temperature,  in  Fig.  10,  yields 

E j -E  j =*  -1.2S  ^ .26  kcal/mole 

5^-3^  > -3.4  + .9  eu  . (A) 

3.  CorracCions 

The  extanc  Co  which  che  ratios  ?,/?.,  cakar.  ac  these  low  temperatures 

L 0 

were  distorted  during  the  quench  can  be  ascimated,  and  corrected  for,  by 
caking  advantage  of  che  race  studies  in  Sec.  III.  We  noce  chac  che  racio  of 
che  jump  races  y^  co  2^,  which  is  required  for  such  studies,  is  equal  co  R^. 
This  racio  is  calculated  from  Sq . (16)  using  che  r3v  daca  in  (A).  The  race 
of  diffusion  of  free  diners,  y,  as  well  as  che  jump  race  3,  from  configura- 
tion 1 co  0,  3ra  sec  equal  co  che  experimental  values  for  rhenium  on 
W(2!L).  Transitions  from  bound  co  dissociated  pairs  are  assumed  co 

occur  ac  che  race  y,  = 2a  exo  - 3500/R,T,  where  R,  is  che  23s  constant. 

Using  chese  race  constants  and  the  methods  outlined  in  3ec.  Ill,  we  can  now 
estimate  the  changes  in  P^/P^  during  che  quench  preceding  field  ion 
microscopy.  The  measured  values  of  ?./?^  are  correcced  for  chese  changes, 
and  raploccad  co  give  a new  estimate  of  che  entropy  and  energy  differences 


between  scaces  0 and  1. 


Another 


le  of  correstions  is  then  initiated, 


using  che  new  parameters  to  evaluate  3 /S^,  and  this  process  is  repeated 
until  further  correction  does  not  significantly  affect  the  energy 
difference  estimated  for  rhenium  dimers  in  state  0 and  1.  As  already 
noted,  only  che  ratio  of  the  rate  constants  3q  to  3^,  not  their  absolute 
value,  affects  the  extant  of  distortion  during  a quench;  this  makes  the 

Q 

iteration  simple.  For  ? /P^,  corrections  even  at  327  K amount  to  only 
107.  and  one  iteration  brings  the  internal  energy  difference  to  within  17. 
of  che  final  value. 

Measurements  of  ? /P^,  corrected  in  this  way,  are  plotted  in  Fig.  10. 
The  thermodynamic  parameters  for  rhenium  dimers  derived  from  the  corrected 


measurements  are 


E^-Eq  = -1.40  + .26  kcal/mole 
S^-Sq  = -3.3  + .9  eu 


The  difference  in  energy  between  state  1 and  0 found  here  is  considerably 
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larger  than  previously  reported  for  the  (211)  plane  of  tungsten.  That  is 
expected;  earlier  results  were  based  on  measurements  extending  to  higher 
temperatures,  at  which  redistribution  during  quenching  dominated.  If  only 
data  taken  over  lower  temperatures  are  compared , then  the  earlier  measure- 
ments agree  with  the  present  to  well  within  one  standard  deviation. 

V.  EQUILIBRIUM  BETWEEN  30UND  AND  DISSOCLAISD  RHENIUM  DIMERS 
With  the  thermodynamics  of  states  0 and  1 in  hand,  the  ratio 
R_  = (?1/?Q)r  can  be  extrapolated  to  temperatures  at  which  dimers  are 
dissociated  to  an  appreciable  extent.  Observations  of  the  frequency  with 
which  dimers  are  dissociated,  that  is  of  M^,  combined  with  relation  (14), 


energy  of  dimers.  One  probLem  sCilL  remains  - Che  racio  P^/P^  as  given  by 
Eq.  (4)  depends  upon  the  Length  L of  the  channel  on  which  measurements 
were  made.  A scheme  of  reducing  the  various  experiments  to  a standard 
length  is  presented  in  Sec.  A.  This  makes  it  possible  to 
analyze  observations  on  the  dissociation  of  rhenium  dimers,  which  are 
presented  in  B.l.  Finally,  corrections  for  the  distortions, which  occur 
on  quenching  to  the  temperature  at  which  field  ion  microscopy  is  done, are 
worked  out  in  Sec.  B.2,  to  yield  the  thermodynamic  parameters  for  Re?  on 
W(211) . 

A.  Normalization  of  Data  for  Channel  Length 
Measurements  of  are  made  on  channels  of  different  length,  and  it  is 
first  necessary  to  reduce  them  to  a common  channel  length  L^.  For  final 
analysis  of  the  data,  we  require  the  ratio  It  is  not  desirable  to 

calculate  this  ratio  for  each  experiment  and  then  to  combine  the  individual 
ratios,  appropriately  adjusted  for  differences  in  L and  weighted, into  one 
final  average.  Instead,  it  is  preferable  to  first  combine  values  of 
and  >1^  from  different  experiments  into  one  grand  total,  prior  to  taking 
the  ratio 

To  accomplish  this  we  take  advantage  of  Eq.  (4)  to  write 
C(L-2)/2](M1/Md)l  = exp  - WL/kT  ; 

here  the  subscript  L explicitly  indicates  the  length  of  the  channel  in  a 

particular  set  of  measurements.  The  left  side  is  now  independent  of 

channel  length,  and  it  is  clear  that  D , the  ratio  M />L  referred 

L0  1 ° 

to  a channel  of  standard  length  Lq,  can  be  deduced  from  experiments 


on  channels  of  different  lengths  L according  to 


31 


■ 


<WLn  • :a-2),(V2)-<VMD)L = °Ln  • (17) 

0 0 

As  already  indicated  La  Sec.  Ill,  Che  racio  (>1^/>L)T  is  noc  obcaiaed  from 
measuremeacs  of  >1^  and  ac  che  same  temperaCure . and  are  derived 
separately,  directly  from  dimer  distributions  ac  high  temperatures,  M_, 
by  extrapolation  from  data  at  lower  temperatures.  However,  their  ratio  is 
readily  corrected  for  the  plane  size. 

Normalization  is  accomplished  by  an  approximation.  We  write  Eq.  (14) 


in  tne  rorm 


(M  ) - M/(l+L/5_  +1/D.  ) 

1 L0  ^0  “0 


The  racio  R,  is  available  by  extrapolation  from  low  temoeratures , in 

L0 

accord  with  Eq.  (16);  a value  of  D is  derived  from  data  on  a channel  of 

L0 

length  L using  Sqs.  (14)  and  (17).  For  a channel  of  standard  length  t , 
the  number  of  dimers  which  are  dissociated  or  in  state  1 can  therefore  be 
obtained  from  che  measurements  using  the  relations 


0 \ 1 


L 

• 0 \ 


I'M  ) = Vf  / 1 4.  — I — 4.  I — -S_ 

^ ' 1 ' - -l/  o 1 1 V _ 1 n 

L Ln  - 'Lfl  L/  *■»  M - - Ur  - 


(Vl  -Ccv2>/<i-*>J<V,/Dt.  • 

o 0 

The  values  (M, ).  . and  (M_ ).  . so  derived  from  the  data  ii 

1 L0,J  73  \),J 

experiment  are  summed  to  finally  give  the  desired  racio 


<VVl  * E (mi}l  i/e  (Vl  i 


Analysis  of  che  sea  rise  Leal  uncerta  inties  in  this  estimate  is  complicated , 
but  follows  the  procedures  usual  in  the  propagation  of  errors,  with  R^, 
the  ratio  ?,/?  on  an  infinite  channel,  and  (M^)^  as  ir.deoendent 


variables . 


3.  Dissociation  Equilibrium 


1.  Observations 


The  equilibrium  between  bound  and  dissociated  rhenium  dimers  on  che 

(211)  plane  of  tungsten  has  been  studied  at  temperatures  ranging  from  351° 

to  412°:<.  At  temperatures  lower  chan  that,  equilibration  is  inconveniently 

slow;  at  higher  temperatures,  the  race  ac  which  adatoms  are  lost  from  the 

0 

plane  is  excessive.  In  the  vicinity  ot  -+00  K,  however,  *»  20*1  of  the  dimers 
are  dissociated,  and  the  equilibrium  between  bound  and  dissociated  pairs  is 
easy  to  observe.  After  depositing  a pair  of  rhenium  adatoms  in  adjacent 
rows,  they  are  equilibrated  at  che  temperature  of  interest  for  a period 
ranging  from  2 co  10  sec.  This  interval  is  chosen  so  that  che  adatoms  in 
a dissociated  dimer  travel  a distance  comparable  to  the  channel  length 
and  can  therefore  sense  each  other's  presence.  The  state  of  the  dimer 
is  then  examined  after  cooling  to  *=  20  K.  This  sequence  :s  repeated, 
until  one  of  the  atoms  is  lost  ac  che  edge  of  the  plane,  or  uncil  a total 
of  roughly  30  observations  has  been  completed.  Thereafter,  the  surface 
is  field  evaporated  for  general  cleanliness,  rhenium  atoms  are  again 
evaporated  onto  che  plane,  and  another  measuring  cycle  is  begun.  Observations 
of  dissociating  dimers  are  show,  in  Fig.  11;  the  overall  results  are 
summarized  in  Table  II. 


(a)  (b)  (c)  (d)  (e) 


(k)  (I)  (m)  (n)  (o) 
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Fig.  11.  Dissociation  sequence  for  rhenium  dimers  on  W(211).  Observations 
after  equilibration  for  3 sec  at  392°K  are  in  chronological  order. 
Labelling  of  atomic  positions  begins  at  left;  channel  length 
L ■ 27.  (a)  Position  of  adatom  in  top  row  - 29,  position  of 

adatom  in  bottom  row  - 19,  configuration  - 10,  (b)  15,  16,  1 
(c)  23,  21,  1 (d)  11,  8,  3 (e)  4,  5,  1 (f)  17,  16,  1 (g)  6,  5,  1 
(h)  4,  11,  7 (i)  4,  4,  0 (j)  12,  11,  1 (k)  15,  15,  0 (1)  2,  2,  0 
(m)  8,8,  0 (n)  27,  26,  1 (o)  14,  25,  11.  In  quantitative 
analysis  of  the  data,  (1)  not  used  as  too  close  to  edge. 
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Measurements  were  made  on  channels  ranging  in  size  from  L = 17  to  32. 

The  data  are  therefore  normalized  to  a standard  channel  length,  as  outlined 

in  the  previous  section;  we  adopt  the  value  LQ  = 25,  which  is  the 

mean  for  our  experiments.  A semi- logar i thmic  plot  of  (P  /P  ) , derived 

' 1 D Lq 

from  Eq.  (21),  versus  1/T  now  yields  the  thermodynamics  of  dissociation, 
in  accord  with  Eq.  (4).  As  stressed  earlier,  the  ratio  P^/P^  not  only 
involves  measurements  of  dissociation,  but  also  of  the  relative  occurrence 
of  states  1 and  0,  obtained  at  lower  temperatures.  If  we  combine  the  raw 
data  on  dissociation,  in  Table  II,  with  values  of  M^/Mq  extrapolated  from 
data  set  (B),  then  we  obtain  the  raw  dissociation  plot  shown  in  Fig.  12. 
This  yields 


E^  = -3.3  + 1.1  kcal/mole 
S1=-1.6+2.9eu 


(C) 


Only  the  statistical  uncertainties  arising  from  measurements  of  M^/Mq  and 
of  Mq , the  number  of  dissociated  dimers,  have  been  included  in  the  error 
bars.  There  is  some  uncertainty  as  well  arising  from  possible  errors  in 
the  lengths  L of  the  channels.  This  uncertainty  is  small,  however, 
amounting  to  only  1%  of  the  standard  error  in  the  internal  energy.  The 
particular  choice  of  Lq  to  which  the  data  are  normalized  also  has  little 
effect;  a change  of  Lq  by  5 alters  the  estimate  of  the  energy  by  only  1". 

2.  Corrections 

The  raw  values  for  the  thermodynamics  of  rhenium  dimers  in  state  1 
are  comparable  to  those  of  the  model  for  which  rate  estimates  were  made 
in  Section  III.B.  For  that  system, redistribution  during  the  time 
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necessary  for  the  sample  to  cool  from  400°X  to  the  temperature  of  observation 

lowered  the  number  of  dissociated  pairs  by  *»  5”,  quite  a significant  amount. 

Similar  effects  must  be  operating  to  distort  the  dissociation  data  for 

rhenium  dimers.  The  number  of  dissociated  pairs  has  there  been  corrected 

bv  the  iterative  method  used  on  the  ratio  ?,/?,  in  Sec.  17.  3.  As  the  first 

1 0 

steo,  the  ratio  2,/3.  is  adjusted  to  t.ne  values  of  and  3,-3,  in  set  (31 . 

• 3 1 10  10 

The  energy  E,  of  the  model  dimers  ^ left  at  -3.5  kcal/aole,  and  the  rate 


changes  in  ? on  quenching  from  the  actual  equilibration  temperatures  are 
made,  and  the  experimental  data  are  corrected  to  yield  an  improved  estimate 
of  the  thermodynamics  of  dissociation.  The  rate  2,,  which  corresponds  to 
the  transition  from  bound  to  free  dimers,  is  adjusted  to  fit  the  new  values 
of  E,  and  5, , and  another  cycle  of  correction  is  initiated.  A total  of 
three  cycles  is  enough  to  converge  on  the  corrected  values  plotted  in 


Wtlicn  i LV0 


S.  * -3.7  + 1.1  kcal/mo le 


5,  = 2.3  + 2.9  eu 


we  consider  these  to  be  the  preferred  estimates  for  the  thermodynamics 


of  rhenium  pairs  in  the  1 configuration  on  the  w (211).  Two  things  should 
be  noted.  The  parameters  of  state  0 and  1 do  r.ot  sensitively  affect  our 
estimates  for  the  thermodynamics  of  dissociation;  s 107,  change  in  E.-Ej 
causes  less  than  i 1"  change  in  5..  Also,  the  various  corrections  for 
distortion  during  quenching  increase  the  internal  energy  by  13*  at  most; 
this  is  a small  effect  tor.oarad  to  the  statistical  uncer taintv . 


dicers.  Bottom  curve  for  raw  dissociation  data  Cy-axis  ac  right) 
yields  thermodynamic  parameters  in  set  (C).  Top  curve  ,y-axis 
ac  left)  shows  dissociation  data  corrected  tor  distortion 
during  quench;  this  yields  dac3  in  (D).  Error  bars  indicate 
standard  deviation  - due  to  statistical  errors  in  observations; 
curves  obtained  by  least-scuaras  fit  to  points  weighted  bv  v~'r-. 
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SUMMARY 

The  measurements  reported  here  'ar^/e  to  define  the  energetic}  of 
rhenium  dimers  on  the  (211)  plane  Jt  tungsten.  Only  one  i-;umptu  ' enters 
the  analysis  of  dissociation  - that  dimers  in  ccnf . guration  ' r ' er 

can  be  treated  as  completely  dissociated.  The  nualer  of  finer; 
in  this  configuration  is  small,  however.  An  approximate  analy*:*  r 
dissociation  is  therefore  possible,  o sed  on  the  alternati.  tisumpt  1 
that  only  dimers  in  configurations  h igher  than  : ne  c . ,ssc  i.i  ted . Th  s 
changes  the  dissociation  energy  by  less  than  10”,  in  imouut  small  n the 
scale  of  the  statistical  scatter.  The  specific  assumptions  about  the 
range  of  atomic  interactions  thus  do  not  significantly  affect  the  results 
on  dissociation 

The  observations  of  the  relative  occurrence  of  dimers  in  states  1 
and  0 reveal  an  unexpected  effect  - the  entropy  of  the  former  is 
considerably  lower.  To  interpret  this  properly,  it  would  be  desirable 
to  compare  the  entropy  of  configuration  0 and  1 with  that  of  dissociated 
diners.  Unfortunately,  the  scatter  in  the  dissociation  results  is  still 
too  large  for  this.  The  internal  energies  of  dimers  are  in  better  shape. 
Compared  with  dissociated  dimers,  we  have 

E^  = -3.7  + 1.1  kcal/mole 
Eq  = -2.3  + 1.2  kcal/mole 

The  cohesive  energy  of  rhenium  pairs  on  the  (211)  thus  is  small, 
on  the  order  of  magnitude  expected  of  van  der  Waals  interactions.  What 
is  surprising  is  the  difference  between  dimers  in  the  0 and  1 configuration. 


aaaaBaaHH| 


38 


The  energy  of  the  0 state  is  significantly  higher  than  of  the  1;  in  the 

9 9 

termer,  the  Interatomic  separation  is  4.48  A compared  to  5.25  A in  the 

l state.  At  cLose  distances,  repulsive  interactions  can  be  expected  due 

to  dipole-dipole  interactions  between  adatoms  polarized  by  the  substrate. 

For  adatoms  with  a net  dipole  of  1 D (a  value  characteristic  of  tungsten 
13 

adatoms  on  the  (110)  plane  of  tungsten),  such  repulsion  increases  the 

•» 

energy  by  only  — 160  cal/mole  at  4.5  A.  This  is  not  of  the  right  magnitude 
to  account  for  the  observed  energy  difference  between  dimers  in  configuration 


0 and  1.  Moreover,  it  must  be  emphasized  that  other  dimers  , for  which 


the  energetics  on  W(211)  have  not  yet  been  quantitatively  examined,  show 
a strong  preference  for  the  0 state  as  compared  to  the  1.  Such  a pronounced 
chemical  specificity  is  not  expected  for  dipole-dipole  repulsions.  Indirect 
interactions  must  therefore  be  responsible  for  these  distinctive  energy 
differences.  Additional  information  about  these  interactions  should  be 
available  soon  from  kinetic  measurements  underway  in  this  laboratory 
on  rhenium  clusters.  We  defer  further  discussion  till  then. 

The  primary  purpose  of  the  present  study  has  been  to  demonstrate 
that  thermodynamic  data  can  be  derived  by  direct  observation  of  metal 
clusters  at  a surface.  With  this  capability  established,  it  would  now 
be  desirable  to  have  quantitative  information,  similar  to  that  presented 
here  for  rhenium  on  W(211),  concerning  the  behavior  of  chemically  different 
adatoms.  This  should  make  it  possible  to  discern  trends  in  atomic  forces 


depending  upon  the  electronic  structure  of  the  atoms  interacting  at  a 


surtace . 
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TABLE  II.  DISSOCIATION  EQUILIBRIUM  OF  RHENIUM  DIMERS  ON  W(21I) 
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22 
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27 

59 
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7 
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25 

3 

12 

50 

5 

26 

3 

3 “J 

Temperature  decay  parameter  = 3.0  x 10  (deg  sec) 

-4  -1 

measurements,  b » 5.75  x 10  (deg  sec) 


For  all  other 
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The  length  of  Che  support  is  a compromise,  influenced  by  two  competing 
effects:  for  a long  loop,  the  temperature  distribution  is  more 
uniform  at  the  sample,  allowing  more  precise  determinations  of  the 
equilibration  temperature;  however,  cooling  is  more  rapid  for 
shorter  loops. 

This  follows  from  Eq.  (Al)  of  reference  5,  with  rate  constants 
conforming  to  the  conventions  in  Fig.  4 of  the  present  paper. 
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FIGURE  CAPTIONS 

Fig.  1.  Hard-sphere  model  of  dimer  or  (211)  plane  of  tungsten.  Dimer 
movement  is  confined  to  [ill]  channels  of  close-packed  atoms. 

a 

Atomic  spacing  l along  channels-2. 74  A;  distance  between  neighboring 

o 

channels  - 4.48  A. 

Fig.  2.  Configurations  of  bound  and  dissociated  rhenium  dimers  on  W(211). 

Hard-sphere  model  on  left,  field  ion  image  on  right.  Dimer 

configurations,  from  top  to  bottom:  0,  1,  2,  and  4;  the  top  two 

are  bound,  the  bottom  two  free.  Interatomic  separation  of  atoms 

in  dimer  established  by  techniques  in  reference  4. 

Fig.  3.  Cooling  curves  for  sample  support  after  equilibration  at  Tq  = 330 
o 

and  402  K.  Points  indicate  temperatures  obtained  from  resistance 
measurements,  solid  curves  are  least-squares  fit  to  1/T  = 1/T  + bt. 

Shading  shows  temperatures  at  which  dimer  distributions  are  frozen 


Fig.  4.  Rate  constants  for  jumps  of  dimers  in  one-dimensional  diffusion. 

Center  of  mass  positions  are  indicated  by  crosses  and  labelled  with 
italic  numbers.  Arrows  show  atomic  jumps  allowed  between 
nearest  neighbor  positions.  Jump  rates  to  higher  configuration 
are  given  by  a,  to  lower  configuration  by  3;  subscript  labels 
the  starting  configuration. 

Fig.  5.  Schematic  potential  diagram  for  model  system.  Dimers  in  configura- 
tion 2 and  higher  are  completely  dissociated;  their  potential 
energy  is  taken  as  the  zero  of  reference.  Rate  constants  are 

assigned  the  following  values:  an/9  = 0.4  exp  1234. 6/RT; 

U 1 G 

= 3 x a exp-3500/RGT,  3^  = 2.5  x lO^exp-lS , 300/R,T  sec"1-, 

12  -1 

a = 2.95  x 10  exp-19, 800/RT  sec  ; the  last  two  rates  are 

G 


4 

experimental  quantities  for  rhenium  on  W(211). 

R * 1.987  cal  (deg  mole)  \ 

G 
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Fig.  6.  Time  evolution  of  dimer  probabilities  for  model  system  at 

-*00°K.  To  bring  probabilities  to  a common  scale,  arbitrary 

constants  are  subtracted  from  ?,  and  P P,  is  multiplied  bv  a 

0 12 

factor  of  10.  P = 1 - (P.+P.+P.). 

s 0 12 

Fig.  7.  Representation  of  coding  curve  for  sample  support  by  isothermal 
intervals  separated  by  step  function  drop  of  A degrees.  For 
quantitative  estimates  of  distortion  during  cooling,  A ~ 2-4°K. 

Fig.  3.  Change  in  probabilities  P^/P^,  P^/PD>  and  PQ  during  quench  of  model 
dimers  equilibrated  at  T^.  Temperature  decay  parameter 
b « 5.75  X 10  * (deg  sec)  ^ . 

Fig.  9.  Time  evolution  of  P^/?^  for  model  dimers  quenched  after  low 

temperature  equilibration.  Temperature  scale  appropriate  only 
to  T^  ■ 340°K.  All  values  for  b.  *>  5.75  x 10  4 (deg  sec)  ^ . 

Fig.  10.  Temperature  dependence  for  ratio  R of  rhenium  dimers  in  state  1, 
compared  to  state  0, on  W(211).  Raw  data  on  bottom  (ordinate  at 
left)  yields  thermodynamic  data  in  set  (A).  Values  on  cop 
corrected  for  distortion  during  quench  (ordinate  at  right),  give 
data  set  (S).  Error  bars  indicate  standard  deviation  arising 
from  statistical  effects;  curves  obtained  by  least-squares  fit 
to  data  points  weighted  by  R*Yvar  R. 

Fig.  11.  Dissociation  sequence  for  rhenium  dimers  on  W(211).  Observations 
after  equilibration  for  3 sec  at  392°K  are  in  chronological 
order.  Labelling  of  atomic  positions  begins  at  left;  channel 
length  L a 27.  (a)  Position  of  adatom  in  top  row  - 29,  position 

of  adatom  in  bottom  row  * 19,  configuration  * 10,  (bt  15,  16,  1 
^ c 0 23,  21,  1 (d)  11,  3,  3 (a)  4,  5,  1 (ft  17,  16,  1 (g)  6,  5,  1 
(h)  4,  11,  7 (i)  4,  4,  0 (j)  12,  11,  1 (k)  15,  15,  0 (It  2,  2,  0 
(m)  3,  3,  0 (n)  27,  26,  1 (ot  14,  25,  11.  In  quantitative 
analysis  of  the  data,  (,1)  not  used  as  too  close  to  edge. 


yields  thermodynamic  parameters  in  set  (C).  Top  curve  (y-axis 

at  left)  shows  dissociation  data  corrected  for  distortion 

during  quench;  this  yields  data  in  (D).  Error  bars  indicate 

standard  deviation-  due  to  statistical  errors  in  observations; 

2 2 

curves  obtained  by  least-squares  fit  to  points  weighted  by  y /a  . 
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